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Abstract We study two different physical scenarios of thermonuclear reactions in stellar plasmas proceeding 
through a narrow resonance at low energy or through the low energy wing of a wide resonance at high 
energy. Correspondingly, we derive two approximate analytical formulae in order to calculate thermonuclear 
resonant reaction rates inside very coupled and non ideal astrophysical plasmas in which non-extensive 
effects are likely to arise. Our results are presented as simple first order corrective factors that generalize 
the well known classical rates obtained in the framework of Maxwell-Boltzmann statistical mechanics. As 
a possible application of our results, we calculate the dependence of the total corrective factor with respect 
to the energy at which the resonance is located, in an extremely dense and non ideal carbon plasma. 

PACS. 24.30.-v Resonance reactions - 26.50.+X Nuclear physics aspects of novae, supernovae, and other 
explosive environments - 05.90.+m Other topics in statistical physics, thermodynamics, and nonlinear 
dynamical systems 



1 Introduction 

Cussons, Langanke and Liolios pQ proposed, on the basis 
of experimental measurements at energy E ~ 2.4 MeV, 
that the resonant behavior of the stellar 12 C + 12 C fu- 
sion cross section could continue down to the astrophysical 
energy range, thus leading rise, beside the usual Debye- 
Hiickel screening (whose corrective factor is /g), to 
further plasma resonant screening effect (conventionally 
described by a suitable /rs factor). While / s > 1 en- 
hances the reaction rate, it has been pointed out that 
Jrs < 1) i- e - the resonant screening effect is likely to 
reduce the rate. The reduction of the resonant rate due 
to resonant screening correction amounts to 11 orders of 
magnitude at the resonant energy of 400 keV, influencing 
the carbon ignition density in white dwarfs. Itoh et al. 3 
have shown that, using an effective screening potential ob- 
tained by one-component-plasma (OCP) Monte Carlo ex- 
periments, the overall effect does indeed strongly enhance 
the carbon-carbon reaction rate by a considerable amount 
(i.e. /t = fs ■ /rs ^ 1) because of the global screening 
domination that amounts to an enhancement of the rate 
by 12 orders of magnitude, with important implications 
for hydrostatic burning in carbon white dwarfs. Given the 
12 q _|_ 12 q reac tion, the current hypothesis is that the 
entrance channel width is much smaller than the total 
resonance width, the latter being much smaller than the 
resonance energy. The same picture could possibly apply 
to other fusion reactions between medium-weighted nuclei 
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(e.g. to the 16 + 16 reaction, which could be an active 
burning stage in some white dwarfs 0]). 

The previous discussion refers to extremely dense stel- 
lar plasmas, characterized by a temperature T ~ 10 8 K, 
a mass density p ~ 10 9 g • cm -3 and a plasma parameter 
r < 178. In these physical conditions, we expect that 
non-extensive effects could also arise. We briefly recall 
that such plasmas show deviations from the several as- 
sumptions that are the basis of Maxwell-Boltzmann dis- 
tribution. Long range many-body nuclear correlations and 
memory effects, among others, can be sufficient to justify 
the use of a distribution function which slightly deviates 
from the standard Maxwell-Boltzmann one [5] (see also, 
for example, ref. |S] for a discussion on physical conditions 
in which non-extensivity needs to be taken into account). 
Our aim is to derive a simple first order formula in order 
to express non-extensive corrections for reactions proceed- 
ing through narrow resonances: we will consider this case 
in subsection 

In addition to the above resonant reactions in white 
dwarfs, let us mention few resonant reactions occurring 
in the stellar interior (like the Sun interior), where no 
carbon burning is active. It is well known that, along 
with the reactions of the proton-proton chain, many re- 
actions of the CNO cycle are presently believed to be 
non-resonant at the relevant stellar energies around the 
maximum Eq of the Gamow peak (whose order of mag- 
nitude is Eq ~ lOkeV inside the Sun's core). However, 
the 17 0(p, a) 14 N reaction, which belongs to the CNO-II 
subcycle, shows a very narrow resonance ( 18 F being the 
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compound nucleus) located at an energy between 65 keV 
and 75 keV, with a proton partial width r p ~ 22 neV; this 
might be taken into account when solar model calculations 
are carried out (see and references therein). Moreover, 
one could actually point out that many non-resonant reac- 
tions in the CNO cycle proceed indeed through the low en- 
ergy wing of a wide resonance located at energy En 3> Eq , 
as was already observed by Clayton |S] (among others, we 
recall also the 14 N(p, 7) 15 reaction; its S factor has been 
very recently measured by the LUNA group jHJ). In this 
case, we could adopt a very reliable resonant formalism 
that allows us to evaluate the astrophysical factor of the 
reaction at the Gamow peak energy, S(E ), by relying ex- 
clusively on the experimental data of the resonance at the 
relatively high energy En- 

For the sake of completing this discussion, let us con- 
sider the 12 C(p, 7) 13 N reaction (CNO-I subcycle). At its 
Gamow energy in the Sun, Eq = 24.68 keV, this reac- 
tion seems to be non-resonant |10) . but it shows a wide 
resonance at energy En — 424 keV Eq, with a total 
width Pi ~ 40 keV, and an electromagnetic channel width 
r i ~ 0.77 eV. It is an easy matter to express the astro- 
physical factor S(Eq) by means of En, It and I 1 -, only, 
through a Breit-Wigner approximation formula 

As it has been pointed out (see, for example, refs. [T2l 
a very slight non-extensive deformation could also 
arise in the Sun interior: we will derive a first order non- 
extensive corrective factor for nuclear reactions proceed- 
ing through the low energy wing of a wide resonance in 
subsection 13. 21 

From the previous discussion it follows that it is useful 
to develop an analytical theory dealing with reactions that 
proceed through narrow or wide resonances; this work was 
already achieved in ref. [BJ. In this paper we extend the 
known classical results to the case in which non-extensive 
corrections arise, and we present a possible application to 
the 12 C + 12 C reaction in a white dwarf's plasma. 



2 Classical resonant reaction rates 

In this section, we briefly review some of the known re- 
sults about the analytical reaction rate calculations in as- 
trophysical plasmas, that were obtained in the framework 
of the classical Maxwell-Boltzmann statistics (MB) [8"lll0| . 

Let us consider a thermonuclear reaction between two 
nuclei i and j. We define the classical reaction rate, r^ B , 



by 



MB 



NjNj 



(Vijdij) MB 



1 



J v Jo 



+oc 



b M B{E)v lJ (E)a ij (E)dE , (1) 



where Ni, Nj are the particle densities, E is the relative 
energy (in the center of mass frame of reference) at which 
the reaction occurs, Vij is the relative velocity between 
two fusing nuclei, 07,- is the reaction cross section, 4>mb is 
the Maxwell-Boltzmann energy distribution function, and 
Sij — 1 if nuclides i and j are identical, 5ij — otherwise. 



As far as resonant reactions are concerned, we can dis- 
tinguish two different cases of physical interest. Let be 
En the resonance energy, Eq the Gamow energy and It 
the total resonance width. If the following conditions are 
satisfied, 

/ E R 1=3 E ,„-. 

I A- « E R ' W 
the reaction proceeds through a narrow resonance at low 
energy (this case will be labelled as "r"). The 12 C + 12 C 
reaction occurring inside a white dwarf, beside the hy- 
potheses already discussed in sect.^ belongs to this first 
scenario. 

On the contrary, if the following inequalities are satis- 
fied, 

j E R > Eq or E-r » Eq 
I It > Eq 



(3) 



the resonance is said to be a wide resonance (this case will 
be labelled as "R"). The 12 C(p, 7) 13 N reaction inside the 
Sun belongs to this second scenario. 



2.1 Narrow resonances at low energy (MB,r) 

If the conditions in eq.J2J hold, the approximate reaction 
rate reads [S] 

r MB,r = (Q T )3/2 N i N j ^ij^ 



1 + S tj (^k B T) 3 / 2 



rin(EYi)r ol 



■ exp 



En \ 



(4) 



where r^En) is the entrance channel width, r ou t is the 
exit channel width, It ~ An + E out is the total resonance 
width, keT is the plasma thermal energy and is the 
reduced mass of the two-body system i + j. In eq.Q, the 
quantum factor u>tj is defined by |14| 

2J C + 1 



(2J i + l)(2J 7 + l) ' 



where Jc is the total quantum angular momentum of the 
compound nucleus and Ji , Jj are the spin numbers of the 
interacting nuclei i and j. 



2.2 Wide resonances (MB,R) 

If the physical conditions in cq.(|3J hold, the approximate 
reaction rate now reads |S] 



mb,r. _ 2 3 / 2 ^ NjNj u; t3 h 2 E„ 



Ql/2 1 i A . . 3/2, m 



r in r ( 



in-* out 



(Eq - En) 2 + (r T /2) 2 



1 i^'- (5) 



Equation (J5J has been obtained by writing the astro- 
physical factor Sij as 



Sij(E) 



TT h 2 



2 ^j (E - E R ) 2 + (r T /2) 2 ' 



(G) 
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through a well known Breit-Wigner approximation for- 
mula, which is supposed to hold if the energy E is close 
to the resonance energy Er. 



3 Non-extensive resonant reaction rates 

In refs. |5llHlll3| it has been shown that a coherent model 
describing a given stellar core might deal with slightly cou- 
pled plasmas, non-Markovian random walks, many-body 
collisions and memory effects. From this point of view, the 
classical theory founded on the Maxwell-Boltzmann statis- 
tical mechanics should be considered a first order approx- 
imation (a very good one indeed inside the Sun's core). 
A further approximation, to which we are referring from 
now on, lies on the more general picture of non-extensive 
statistical mechanics |15Hlfi) . 

Here we limit ourselves to briefly outline, by means 
of three different approaches, how non-standard statis- 
tics can arise. In Fokker-Planck equation context, we can 
introduce corrections to the friction and diffusion coef- 
ficients, considering their expressions to the next order 
in the velocity variable; stationary solution is the Tsal- 
lis non-extensive distribution. In a plasma, each particle 
is affected by the total electric field distribution of the 
other charges; the total microfields have a relatively small 
random component, generally show long time correlations 
and generate anomalous diffusion. All these effects imply a 
deviation from the Maxwell-Boltzmann distribution whose 
entity depends on the plasma parameter and on an ion-ion 
correlation parameter. Correlations among the collective 
modes can lead to a long time asymptotic behavior of the 
velocity correlations of the ions and to anomalous diffu- 
sion, that are related to generalized entropy and generate 
non-Maxwellian probability distributions. 

In this new picture, the reaction rates are defined in the 
same way as in eq.Q, but now the non-extensive energy 
distribution function </>ne (E ) is used, instead of the usual 
Maxwell-Boltzmann distribution 0mb(^)- 

It has been shown |6lll2lll3| that the analytical rela- 
tionship linking </>ne and </>mb can be cast, to first order 
of approximation, in the following fashion 



(E) exp 



-6 



E 
kr/F 



(7) 



provided that \S\ < k B T/E (or \6\ <C k B T/E 0l for com- 
putational purposes). The order of magnitude of the de- 
formation parameter S is 10~ 3 -r 10~ 2 inside the Sun, 
but it can be higher in extremely dense stellar plasmas 
(e.g. in white dwarfs). The 6 parameter is linearly re- 
lated to the q entropic parameter that appears in the non- 
extensive formalism |15) : the relationship between them 
two is 6 = (1 — q)/2. 

We can distinguish two physically relevant cases: if S < 
0, the high energy tail of the </)ne distribution function is 
increased with respect to 4>mb (super-extensivity) , while 
if 5 > 0, the high energy tail is depleted (sub-extensivity) . 

Now we are able to develop the entire non-extensive 
resonant formalism of the thermonuclear reaction rates. 



3.1 Narrow resonances at low energy (NE,r) 

From eqs.Q and (JJJ), the non-extensive rate of a given 
reaction proceeding through a narrow resonance at low 
energy reads 



NE,: 



NiNj 



+ OC 



x v ij (E)a ij (E)exp 



4>mb{E) 

2' 



E 



k R T 



dE 



(8) 



Starting from the hypotheses in eq.J2J we can state 
that, in the energy interval Er — Tt < E < Er + Tt, the 
integral in eq.© is strongly ruled by the reaction cross 
section function <xy(-E) only. Thus we can write, with very 
good approximation, the following result 



NE,r 



r MB ' r [ 1 



exp 



-S 



k^T 



(9) 



If \S\ <C (keT/i^ri) 2 , we can linearize eq.(0, and the 
first order formula is 



NE,j 



MB,: 



[l + Ci(k B T,£ R )<S] , 



(10) 



where 



<*<W.*0-T-(& 



Therefore, our final result in the case of narrow reso- 
nances at low energy is expressed by 



NE.r MB.r 
f . . ' = f . . 



1 + 15 A f Er 



If the condition 



En 



E Q 



k R T k R T 



> 



1. <):>>(; 

4 



(11) 



(12) 



is satisfied, then Ci(kBT, Er) < 0. In the Sun interior, the 
thermal energy is k R T ~ 1.36 keV and Eq is in the energy 
interval 24keV < E < 29keV for the CNO cycle, and 
Eq w 6 keV for the p + p — > d + e + + v e reaction. Thus, 
in the Sun's core, and in many other stellar plasmas of 
interest, eq.lJEJ is always satisfied; then, from ea. flfUl . we 
can actually state that the non-extensive rate is increased 
or diminished with respect to the classically calculated 
rate, whether 5 < or <5 > (as already mentioned in 
sect. El • 



3.2 Wide resonances (NE,R) 

In the case of a generic reaction (no matter if resonant or 
not), the following formula linking r^ E and r^ B holds, 



NE MB 



Sij(E Q ) 



13 S l3 {E ) 



15 , 



3 k B T/ 



(13) 
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where 

1 - 




-2/3' 



and 



E = E 1 + 26 



Eg 
k B T 



-2/3 



,(14) 



(15) 



Equations and (|14|l express the first order cor- 
rection to the classical reaction rate r^ B , due to the de- 
formed distribution function that has been already defined 
in eq.0. Besides, ca. H15(l shows implicitly the relationship 
between the new (deformed) Gamow energy Eq, and the 
classical one Eq. A complete proof of these equations can 
be found in ref. |12| . in which the authors started from an 
ad hoc assumption of the deformed distribution function, 
0ne(B), regardless of any statistical basis. In this paper, 
we adopt the same results, but lying on the ground of 
non-extensive statistical mechanics (as was outlined, for 
example, in ref. 

In ea. lfT^jl . the physical properties of the reaction are 
summarized in the ratio between the astrophysical factor 
at the new Gamow energy, Sij(Eo), and the same factor 
at Eq energy, Sij(Eo). Our aim is to express that ratio as 
a function of Br and Br only. The deformation parameter 
5 is assumed to satisfy the following inequality 



k B T 
E 



Then we can write, as a formal expansion in 5, 



(16) 



(17) 



where the two coefficients S and Si can be calculated as 
follows. From eq.©, we immediately obtain that 



So — Sij(Eo) 



it h 2 



^ln-^out 



2 Mij (B - B R ) 2 + (B T /2) 2 



(18) 



On the contrary, the first order coefficient in ea. (|17J) 

reads 



Si = TTUJij 



(Eq — En)r in r out dB 



My [(B -B R ) 2 +B 2 /4] 2 dS 



(19) 



5=0 



In order to calculate the analytical expression for the 
dEo/dd\s=o derivative, we differentiate both members of 
ea. <|15|) with respect to 6, obtaining 



dE 
dS 




-5/3 



3k B T 



from which, without any approximation, it follows that 

(20) 



dEo 
dS 



s=o 



_____ 

3k R T 



Now we can rewrite ea. l|17[) . using the results of eas. l|18fl 
and l|_D|) . as 

S^(E )^S t3 (E a )(l + ^S 

and then, the Sij(Eo)/Sij(Eo) ratio eventually becomes 
Sij{E ) , , 8 r E a — En 



Sij(Eo) 



1 



El 



3 (£ -£ R ) 2 + r|/4k B T 



(21) 



Placing the previous result of eq.JSTJ into ca. (fTTT|l . it is 
clear that the non-extensive reaction rate, to first order, 
reads 



NE,R MB,R / A \ 



1 



8^ (E Q - E R )E 



4 3 k R T 



E, 







3 (Eq — -Br) 2 + -Tr 2 ,/ 4 k B T 



(22) 



A further approximation, under the hypothesis stated 
in ea. lfT?))l . is 



cxp(-Z\ ii ) « 1 



Eq 
k R T 



and therefore, from ea. H22|) . the final result immediately 
follows 



NE,R M B ,R 

r ■ — : r ■ ■ ' 



4 3 k B T 



f_i.V 



i 8 (Eq - Er)Eq 



Eq 



3 (E - Br) 2 + T 2 /4 k B T 



(23) 



4 Conclusions and discussion 

In this work we have analytically derived two first order 
formulae that can be used to express the non-extensive 
reaction rate as a product of the classical reaction rate 
times a suitable corrective factor for both narrow and wide 
resonances, as shown in eas. ljlH and l|23|) . It should be 
stressed that the previous results are correct only if \S\ 
is very small (in the sense of ea. l|16|) V namely if we are 
dealing with slight deformations of the energy distribution 
function. This is really the most common situation, as 
far as astrophysical plasmas are concerned (in fact \S\ ~ 

icr 3 ^i(r 2 ).' 

Concerning the fusion reactions between two medium- 
weighted nuclei, for example the 12 C + 12 C reaction, our 
non-extensive factor, which now can be formally defined 
as follows, 



NE 



15, _(____ 
4 U B T 



gives rise to further correction beside the screening and the 
potential resonant screening, already investigated in pQ 
and 
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Figure 1. Linear plot of the effective factor F, denned in 
ea. l|24fl . against the resonance energy Er. The dash-dotted 
(upper) line refers to super-extensivity, the dashed (lower) line 
to sub-extensivity, while the solid (middle) line describes the 
classical (MB) result. 

It is important to point out that the non-extensivity 
does not affect the other plasma corrections: therefore, we 
can define an effective factor F as 

F = /ne • fs ■ hs , (24) 

where fs and /rs account for the Debye-Hiickel screening 
and the resonant screening effect respectively. 

We have applied our results to a physical model de- 
scribing a carbon white dwarf's plasma, with a temper- 
ature of T — 8 • 10 s K and a mass density of p = 2 • 
10 9 g • cm -3 (the plasma parameter is, correspondingly, 
r ~ 5.6). Furthermore, we have set a deformation param- 
eter \S\ = 10 -3 , regardless of its sign, and we have kept the 
energy of the possible resonance, Er, as a free parameter. 

In fig. n we plot the total corrective factor defined in 
eq. I|24[) against the energy at which the narrow resonance 
is supposed to be located, in the range between keV and 
2400 keV, considering three very different conditions (de- 
pending on the sign and value of the deformation param- 
eter 5). In our calculations, in order to estimate the func- 
tional dependence of fs ■ /rs with respect to Er, we have 
adopted the fitting formulae derived by Itoh and collab- 
orators |2| : their result, that has been derived through a 
classical treatment, is also shown in fig. H by the line la- 
belled with <5 = 0. From the same figure, it is clear that 
a slight non-extensivity does introduce non-trivial correc- 
tions that become more and more important as the res- 
onance energy rises: if 6 > 0, our total effective factor is 
F ~ 136 at Er = 2.4 MeV while, at the same resonance 
energy, the factor is F ~ 1484 if <5 < 0. Anyway, in both 
cases, lim F = 1 when Er — » keV. It is also noteworthy 
that the effective factor always acts to enhance the res- 
onant reaction rate, no matter if the system is super or 
sub-extensive. In conclusion, all the plasma enhancements 



due to the presence of long range many-body nuclear cor- 
relations and memory effects, that can be described within 
the non-extensive statistics by means of the entropic pa- 
rameter q > 1 (S < 0), are in the direction of still more 
increasing the effective factor F of nuclear rates of hydro- 
static burning and white dwarfs environment. 
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